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Abstract 



Using molecular dynamics simulations we study the dynamics of a water-like TIP5P model of 

o ' 

water in hydrophilic and hydrophobic confinement. We find that in case of extreme nanocofinement 



such that there is only one molecular layer of water between the confinement surface, the dynamics 



a. 

of water remains Arrhenius with a very high activation energy up to high temperatures. In case 
Q of polar (hydrophilic) confinement, The intermediate time scale dynamics of water is drastically 

o ! 

modified presumably due to the transient coupling of dipoles with the effective electric field due to 

the surface charges. Specifically, we find that in the presence of the polar surfaces, the dynamics of 

O , 

monolayer water shows anomalous region - namely the lateral mean square displacement displays a 

CO ' . 

distinct superdiffusive intermediate time scale behavior in addition to ballistic and diffusive regimes. 



We explain these finding by proposing a simple model. Furthermore, we find that confinement and 
the surface polarity changes the vibrational density of states specifically we see the enhancement 
of the low frequency collective modes in confinement compared to bulk water. Finally, we find that 
the length scale of translational-orientational coupling increases with the strength of the polarity 
of the surface. 



I. INTRODUCTION 



Water is the most ubiquitous liquid and hence plays a very important role in different 
phenomena. The list of such phenomena ranges across as varied disciplines as physics, 
chemistry biology, nanofiuidics, geology and atmospheric sciences [IKD]. Water also falls 
into a class of complex liquids which are known to have anomalous behavior compared to 
simple liquids |8|, |9j. The anomalous expansion upon decreasing temperature below 4° C and 
increase of specific heat upon decreasing temperature in the supercooled state are some of 
the examples of anomalies of water. Indeed, the list of the anomalous behaviors is incomplete 
and seems to be ever growing lOj, 111 ] . 

Studies of water confined in carbon nanotubes and between hydrophobic surfaces sug- 
gest that the confining surfaces may induce layering in liquid water in extreme nano con- 
finements 12j, |l3j. Moreover, first order layering transitions are also observed 12|, [lj] in 



hydrophobic confinements. Water can also form wide variety of crystalline phases such as 
monolayer and bilayer structures 12Hl8| in nano confinements. Apart from distinct struc- 
tural changes, water also exhibits dynamic changes depending on the nature of the confining 
substrates. It has been found the dynamics becomes slow near hydrophilic surfaces compared 
to hydrophobic surfaces, presumably due to the binding of polar oxygens and hydrogens to 
the polar groups on the hydrophilic substrates. Moreover, the dynamics in confinement 
may depend on the surface morphology. Specifically, it was found that a liquid in smooth 
confinements diffuses faster than when the confining surfaces are rough 19] . Recent studies 
la . l20M23| suggest that both thermodynamic and dynamic anomalies of water are strongly 
affected when water is confined in small nanopores. It is found that when water is confined 
between hydrophobic surfaces both thermodynamic and dynamic anomalies shift to low 
temperatures and low pressures (l8ll20|. Although the characteristics of hydrogen bond dy- 



namics of water is bulk-like in hydrophobic confinement, the average life-time decreases 24 ] 



leading to a shi 



plane 



18 



20 



21 



t in the regions of anomalous dynamics and thermodynamics in the P-T 



231 ] . Moreover, the long time relaxation of hydrogen bonds depend on the 



effective system dimensionality [24]. 

One of the many important questions that can be asked from the biological science 
perspective is - what role, if at all, water plays in biology ? The answer to such questions is 
probably not very apparent and requires a systematic effort to understand physical behavior 



of water under various conditions such as - water in and around ion channels, protein- 
hydration water, water in different solutions etc. Although there has been a large array of 
works answering some of the questions 3j, \% |25I437||. the behavior of water is far from fully 
understood. 

In this paper, we study the dynamics of a monolayer water confined between both hy- 
drophobic and polar surfaces. In section II, we describe the system and method and discuss 
the results of positive and negative surface polarity on the dynamics of water in section 
III. In section IV, we present the results of the orientational dynamics and translational- 
orientational coupling and finally we conclude with a discussion and summary in section 
V. 



II. SYSTEM AND METHOD 



We perform the molecular dynamics simulation of TIP5P (transferable interaction po- 



tential five points) [38H41J water-like molecules confined between two structured surfaces 
separated by 0.6 nm. The separation is chosen such that there is effectively one layer of 
water between the surfaces (see Fig. [T]). The atoms on the surfaces are arranged in a hexag- 
onal packing with as = 0.339 nm and do not interact with each other. The positions of the 
surface atoms were constrained to their respective mean positions by a harmonic potential 
with spring constant 100 kJ/mol/nm 2 . In the case of hydrophobic surface, the interaction 
between the surface atoms and the oxygen of the water molecule is modeled using 6-12 LJ 
potential U(r) 

U(r) = Ae[(^r-^} (1) 

where r so is the distance between the surface atom and oxygen of water molecule, and 
a = (a s + a ow )/2 and e = ^es^ow- <?nw an d tow are the parameters of the LJ interaction 



between oxygens of water molecules [38L l39| . The hydrophilicity of the surface is achieved 
by adding point charge Qe Coulomb to each surface atom and in that case long-range 
Coulombic interactions between surface atoms and water molecules are added in addition 
to the LJ interaction. We performed simulations of hydrophobic and hydrophilic confined 
water in NVT ensemble with the effective density of water p = 1.00 g/cm 3 . The equations 



of motion are integrated with a time step of 0.001 ps and Berendsen's thermostat is used to 
attain constant temperature. 

III. TRANSLATIONAL DYNAMICS AND VIBRATIONAL DENSITY OF 
STATES 

The translational dynamics in the confined systems can be described by lateral mean- 
square displacement (MSD) (r?,(t)), parallel to surface in the periodic directions. In 2- 
dimension, the diffusion constant D\\ can be calculated by using a modified Einstein relation 
between the long time behavior of the (r?,(t)) and D\\ |42J 

Du = lim -i-— (2) 

11 t^oo At v ; 

We first study the temperature dependence of (rjj(£)) and D\\ in the case when there are 
no charges on the surface. In Fig. [2] (a), we show the MSD (rjj(t)) for different temperatures. 
The MSD shows typical time scales -namely a ballistic regime at very small time scales, a 
cage regime at intermediate time scales, and diffusive regime at long times. In Fig. [2j we 
show the Arrhenius plot of D\\ extracted from the long time linear regime of the MSD. We 
find that the dynamics remains Arrhenius for all the temperature range studied and can be 
fit with D\\ = D\\(0)e~ EA / kBT , where Ea = 14.12 kJ/mol and ks are the activation energy 
and the Boltzmann constant respectively and -Dn(O) is a fitting parameter. An Arrhenius 
behavior of D\\ suggests the the water remains very ordered upto high temperatures in 
extreme hydrophobic confinements. Moreover, the activation energy is comparable to the 
low density liquid (where the local structure resembles the Icelh) of bulk water phase. 

We next study the translational dynamics in the case when surface atoms have charges 
on them. In Figure [3](a) and (b), we show the MSD displacement at T = 300 K as a 
function t for different strength of surface polarity. For the hydrophobic case when Q = 0, 
we find that the dynamics of water in confinement is much slower with a lateral diffusion 
constant D\\ ~ 8.0xl0~ 6 cm 2 /s. Note that the diffusion constant of TIP5P bulk water at the 



same density and temperature is ~ 2.5xl0 _5 cm 2 /s 40|. When the polarity of the surface 
increases, there are namely two changes in the behavior of the MSD - (i) the change in 
the behavior of cage dynamics or the intermediate time scale dynamics and (ii) increase of 



the long time MSD and hence the diffusion constant Dn. In Fig. [3](c) and (d), we show 
D\\ for positive and negative polarity of the surfaces respectively. As we would expect the 
dependence of D\\ on strength of polarity in both the cases are similar. 

Moreover, we find that dynamics in polar confinements is different from dynamics at 
high T in apolar confinement. A comparison of MSD at high T and Q = with MSD for 
T = 300 and Q = 0.06 suggests that while at high T and Q = the MSD changes directly 
from ballistic to diffusive behavior, the MSD shows a distinct difference at the intermediate 
time scales in case of polar confining surfaces. 

The change in the behavior of the intermediate time scale dynamics can be understood by 
the fact that although the time average of external electric field at a point inside the confined 
region is zero, a transient coupling of the fluctuating surface electric field E((x),t) with the 
dipole jl(x,t) of a water molecule at position X will impart a torque fI(x,t)XE(x,t). The 
electric field E(x, t) due to one of the surfaces at a position x is given by 

1 N " - _ - 

*(*■«> = s;5>j£fp (3) 

»=1 ' l| 

where the sum is taken over N q charges on the surface and Xi are the position vectors of 
the point charges on the surface. The dynamics of Xi follows the dynamics of a Brownian 
harmonic oscillator and can be written as: 

-^ = -mu> 2 Xi-tvi + F(t), (4) 

where F(t) is a random force and ( is the frictional attenuation due to water. Note that in our 
system the surface atoms do not interact with each other. Although the set of two equations 
describes the time dependence of surface electric field, it is almost intractable. To simplify 
this, we assume that the electric field inside the confined region is due to a flat surface with 
a fluctuating surface charge density cr(t) and hence E(t) now will be perpendicular to the 
surface and can approximately be written as: 



m = ^ (5) 



where the surface charge density a(t) is 



a{t) = MT) (6) 



where A(t) is the surface area at a given time t. Note that the expression for electric field 
is approximate but we assume that the extent of the surface area considered for a(t) is very 
large compared to the distance of the point at which we are interested to find the electric 

field at. Let's further assume that A{t) = A(l -±1 ), we can write the instantaneous 

fluctuations Sa(t) in a as: 

MO = ^j£W), (7) 

where A is the average area of consideration. Using Eq. [7] we can the time autocorrelation 
of fluctuations in cr(t) as: 

N 2 2 
(Sa(t)Sa(0)) = -g-(6A(t)6A{0)) (8) 

Hence the time autocorrelation in the instantaneous fluctuations SE of E can be written as: 

N 2 2 
(6E(t)6E(0)) = ^(6A(t)6A(0)) (9) 

Now assuming that A = Lx.Ly = (x n — Xi)(y n — y\) = xy, where x%,yi and x n ,y n are the 
coordinates of the I s * and n th point charge along the boundary of the area (note that for 
simplicity we only cosider the in-plane motion of the point charges) . The dynamics of the 
new coordinates x and y are given by Eq. HI (SE(t)SE(0)) can now be written as: 

N 2 Q 2 _ 2 

(6E(t)6E(0)) = ^fL(8x(t)6x(0)) = &j(8x(t)5x(p)), (10) 

where Sx = x — Sx(t) is the instantaneous fluctuation in x(t) with (5x(t)) = O.For simplifi- 
cation in the above equation, we have assumed x = y. Using Eq. |4]for coordinate Sx, the 
time autocorrelation of SE can be written as: 

(SE(t)SE(0)) = E 2 ^^e-^ t (cosco 1 t + -1—sinwit), (11) 

AmuiQ 2mwi 



where Uq = \ — is the characteristic frequency of point charges with spring constant k and 



mass m, and uj\ = \J u)q — ^-j. We can write the time autocorrelation of instantaneous 
fluctuation SErif) in the net electric field can be written as: 



, 2 4 k B T _ c t , , 7 

A muiQ 2mwi 



{SE T (t)SE T (0)) = E 2 ^^^e-^ t (cosco 1 t + — r —sinuit), (12) 



From Eq. [T2J we can see that the fast relaxation of autocorrelation of fluctuations in total 
electric field at a point in the confined region occurs with relaxation time 2m/ (. One would 



expect that over this time scale the net electric field would not completely vanish in the 
confined region due to both surfaces and hence there would be a transient coupling of water 
dipoles with the net electric field. If this energy of dipole coupling is comparable to small 
frequency collective modes, such transient coupling of dipoles with the surface electric field 
will lead to a change in the intermediate time scale dynamics. Assuming a surface with 
uniform charge distribution, a simple estimate shows that Q = 0.10 would be equivalent to 
an effective electric field barE « 10 8 V/m. Assuming the length of the considered area to be 
twice the cut-off length (« 2nm), the standard deviation of the net electric field is ~ 0.22E 
at T = 300 K. In this case, the field-dipole coupling would be large enough to change the low 
frequency collective vibrations. To verify our picture we calculate the vibrational density of 
state D(w) from the Fourier transform of the velocity autocorrelation function C(t) defined 
as 

poo 

D(u) = / C(t)e iujt dt, (13) 

Jo 

where u is the frequency and C(t) is defined as 



N 



where Vi(t) denotes the velocity of the i th -water molecule at time t and angular brackets 
denote the ensemble average. In Fig. |5fa) and (b), we show C(t) for different values of 
positive and negative polarity of the surface respectively. We find that as the magnitude of 
the polarity increases, the anticorrelation in C(t) at t ~ 0.03 ps decreases. In Fig. [5] (c) and 
(d), we show the vibrational density of states obtained using equation Eq. [131 A comparison 
of the bulk TIP5P and confined vibrational density of state shows that the low frequency 
collective modes are enhanced in confinement. Furthermore, the hydrogen bond stretch 
mode (~ 210 cm -1 ) is absent (or not distinguishable) in the extreme confinement case (see 
Fig. [5]). We find that when magnitude of the polarity increases, the peak corresponding 
to the OOO bending mode (~ 50cm _1 ) gets broader and merges into the diffusive modes 
suggesting that highly charged surfaces can modify low frequency collective modes and hence 
may lead to destruction of the local order. 



IV. ORIENTATIONAL DYNAMICS 

We next study the orientational dynamics of water molecules since we expect that the 
transient coupling of fluctuating surface electric field to dipoles of water would provide 
exteral torque which would tend to align the water dipoles. To study the orientational 
dynamics of water, we calculated the self-dipole orientational correlation function Cd(t) 42| 
defined as 

where jli(t) is the dipole of i th water molecule at time t, cos 9i(t) is the angle between the 
dipole vectors of i th -molecule at time t and time 0. 

In Fig. [6] (a) and (b), we show Cd(t) for different positive and negative values of Q. We find 
that, for small magnitude of Q, the orientational correlation function Cd{t) decays to zero 
slowly. To quantify the characteristic time of orientational correlations, we define correlation 
time t as the time at which Cdit) decays by a factor of e. In Fig. O (c) and (d), we show r 
for the positive and negative surface polarity respectively. We find that, as the magnitude 
of Q increases, r decreases and seems to saturate to a constant value for larger magnitudes 
ofQ. 



A significant coupling of translational and orientational motion is found in bulk water 43L 144 1 
and other confined systems 45[. Since we can expect that the presence of the charges of 
the confining surfaces strongly influence the dipole fluctuations of the confined water, we 
would also expect that changing the polarity may lead to a change in the effective length 
scale of the coupling of the translational and orientational dynamics. The characterize the 
the effective length scale of the coupling between translational and orientational motion, we 
define £ 



£ =V /D||*r (16) 

We interpret £ as the effective length scale over which translational and orientational motions 
are coupled. Figure [7] shows £ for different values of Q for both positive and negative Q. 
The value £ increases as the Q is increased and seems to aysmptotically reach a constant 
value. Furthermore, we find that the value of £ is similar for both positive and negative 
surface polarities. 



V. SUMMARY AND DISCUSSIONS 

In summary, using molecular dynamics simulations of TIP5P model of water, we have 
investigated the dynamics of a monolayer water system confined between hydrophobic and 
hydrophilic surfaces of different polarity. We find that dynamics of monolayer water in 
hydrophobic confinement remains Arrhenius upto very large temperatures with a very high 
activation energy E& = 14.12 kJ/mol. Further, we find that the transient coupling of the 
water dipoles with the electric field of the surface gives rise to distinct cage dynamics. 
Specifically at high polarity the slope of the MSD as a function of t at the intermediate 
time scales becomes > 1.0 on a log- log scale, giving rise to three distinct dynamic regions 
-ballistic at very small time scales, super diffusive at intermediate times and a long time 
diffusive behavior. Note that although the effect of transient coupling of water dipoles to 
the external field is effective increase of diffusion of the system, the dynamic changes are very 
different in the two cases. While at high temperatures MSD would only show two distinct 
regimes, water in polar confinements has a different intermediate time scale dynamics and 
three distinct regimes. We explain these finding by proposing a simple model of electric 
field fluctuations in the confined region. We further find that the extreme confinement 
enhances the low frequency collective modes. Moreover, we find that the presence of a 
highly polar surface disrupts the low frequency collective vibrational modes and hence leads 
to less ordered system, suggesting that while extreme hydrophobic confinement increases 
the order, the presence of a highly polar surface makes the water molecule relatively more 
disordered in confinement. Finally, we calculated the effective length scale of coupling of 
translational and orientational motions and find that the coupling length scale increases as 
the polarity of the surface increases. 
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FIG. 1: (Color online) Schematics of monolayer water confined between two surfaces, (a) The top 
view of the confining surfaces, (b) The top view of the confining surface with the water molecules, 
(c) The lateral view of the confined system, (d) Density profile p{z) of the surface atoms and the 
water molecules along the confinement direction. The effective confinement width due to excluded 
volume interaction between the water molecules and the surface molecules is 0.320nm, which is 
about the diameter of a water molecule. 
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FIG. 2: (Color online) (a) (rn(t)) as a function of t for different temperatures when there are no 
charges on the surface, (b) A linear dependence of Du on 1/T on Arrhenius plot shows that the dy- 
namics remains Arrhenius even at high temperatures with an activation energy Ea = 14.12 kJ/mol. 
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FIG. 3: (Color online) (a) Lateral mean square displacement (r?,(t)) as a function of t for different 
values of Q. Increasing polarity changes the intermediate time scale dynamics - namely the slope of 
the cage region increases as the polarity increases. For large values of Q, the cage regime becomes 
superdiffusive. (b) (r?,(t)} for the case of surface with negative polarity. The behavior in the case 
of negative polarity is similar to the case of positive polarity (see Fig. [3] (a)) (c) Lateral diffusion 
constant Dm of water molecules increases with increasing Q and seems to asymptotically saturate 
to a constant value for large Q. (d) Analog of Fig. Etc) for negative surface polarity. 
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FIG. 4: (Color online) (a) Lateral mean square displacement (r?<(t)) for Q = 0.20 showing three 
distinct regimes in the dynamics - small time ballistic motion, intermediate time superdiffusive, and 
long time diffusive behaviors, (b) Comparison of MSD at high T and Q = with MSD for T = 300 
and Q = 0.06 suggest that while at high T and Q = the MSD changes from ballistic to diffusive 
behavior as expected for bulk liquids, the MSD shows a different behavior at the intermediate time 
scales in the case of polar confining surfaces. 
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FIG. 5: (Color online) (a) Velocity autocorrelation function C(t) as a function of time for different 
values of the polarity of the surface. The anticorrelation in C(t) begins to disappear as the polarity 
of the surface increases, (b) Analogue of Fig. 0a) for negative polarity, (c) Vibrational density 
of states as calculated from the fourier transform of C(t). A comparison with bulk water D(cj) 
shows that the hydrogen bond stretch band (~ 210 cm" 1 ) is absent in the confinement. Note the 
merging of diffusive modes into low frequency collective modes when Q increases, (d) Analog of 
Fig. [5^c) for negative polarity. 
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FIG. 6: (Color online) (a) Dipole-dipole orientational correlation function C(t) as a function of t for 
different values of the polarity. C(t) decays faster as the surface polarity is increased, (b) Analog 
of Fig. [6^a) for negative surface polarity, (c) Orientational correlation time r as a function of Q 
suggests that the orientational correlation time increases monotonically. (d) Analog of Fig. [f^c) 
for negative surface polarity. 
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FIG. 7: Translation-orientational coupling length scale £ as a function of atomic charge on the 
surface. The circles and squares are for positive and negative polarity respectively. £ increases as 
the the magniture of the charges increases for both polarities of the surface. 



18 



